Hard-Sphere Brownian Motion in Ideal Gas : 
Inter-Particle Correlations, Boltzmann-Grad Limit, 
and Destroying the Myth of Molecular Chaos Propagation 
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The BBGKY hierarchy of equations for a particle interacting with ideal gas is analyzed in terms of 
irreducible many-particle correlations between gas atoms and the particle's motion. The transition 
to the hard-sphere interaction is formulated from viewpoint of the recently discovered exact relations 
connecting the correlations with the particle's probability distribution. Then the Boltzmann-Grad 
limit is considered and shown not to lead to the Bolzmann hierarchy and the molecular chaos, since 
correlations of all orders keep significant. in this limit, merely taking a singular form. 
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I. INTRODUCTION 

In the work [l[ , as well as in its arXive preprints [2] and 
premising works [1, 0] , based on principles and methods 
of rigorous statistical mechanics, - first of all, on the Bo- 
golyubov approach to it [51] , - 1 proved existence of general 
exact relations ("virial relations") connecting (i) proba- 
bility law of random walk of a test "Brownian" particle 
(BP) in a fluid and (ii) irreducible many-particle statis- 
tical correlations between molecules of the fluid and the 
bp's walk (in particular, BP may be merely a marked 
molecule). I emphasized also that these relations make 
it clearly visible that all n-particle correlations always 
are quantitatively significant, even under the low-density 
("Boltzmann-Grad") gas limit. This fact, in turn, implies 
that the Boltzmannian kinetics is incorrect even in this 
limit, and true kinetics of fluids, - including low-density 
gases! - should take into account all the correlations. 

The meaning of so universal correlations was explained 
already in [6] (this article is hardly available on-line, but 
one can read its my own author's translation Q from 
Russian original, or see preprint fsll). It is sufficient to 
notice that a fluid as the whole (and first of all dilute gas!) 
is indifferent to a number of past collisions happened to 
its given particular molecule. Therefore any molecule has 
no definite ( a priori predictable) "time-average" rate of 
collisions 0, In other words, actual rate of its colli- 
sions undergoes fluctuations which are indifferent to time 
averaging. Thus that are scaleless fluctuations with 1/f 
type spectrum fTT'|. The mentioned statistical correla- 
tions directly reflect complicity of particles (via mutual 
collisions) in these fluctuations and, hence, indirectly de- 
scribe their statistics (for detail see ^ or ff] or Q). 

As far as I know, first such statements about molec- 
ular motion were put forward in works by G.Bochkov 
and me [T2| - [l7j as conjectures about origin and prop- 
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erties of 1/f noise accompanying charge transport (i.e. 
Brownian motion of charge carriers) and other transport 
processes (in generahzed sense We demonstrated 

once again that correct ideology leads to useful results 
even at phenomenological level. In particular, in fl^i llBl 
it was shown that fluctuations of rate of collisions (or, 
equivalently, 1/f fluctuations of diffusivity and mobility) 
possess essentially non-Gaussian statistics gravitating to- 
wards power-law probability tails. 

The latter circumstance requires, in view of the well 
known Marcinkiewicz theorem (see e.g. [3), to deal 
with whole infinite chain of cumulants of the fluctuations. 
On the other hand, according to the later works made at 
the microscopic level, firstly i6j| and then its development 
in 19|, n-order cumulant of fluctuations in the rate of 
collisions associates with specifically (n-fl)-particle cor- 
relations (irreducible component of (n + l)-particle dis- 
tribution function) . That is why one should not truncate 
the BBGKY hierarchy of equations! Anyway, neglect of 
three-particle (and thus higher) correlations means rejec- 
tion of the fiuctuations at all (like neglect of two-particle 
correlations rejects any collisions at all [20|). 

The aforesaid was confirmed by exact results of [iHl]. 
Nevertheless, a specially visual analysis of inter-particle 
correlations may be useful. This is just one of purposes 
of the present paper. 

With this purpose it is natural to concentrate on the 
Brownian motion of (molecular-size) particle in ideal gas 
(see [m or some of preprints [12] and also [H, 0]), which 
is most simple "kinetic process" since it produces least 
amount of inter-particle correlations. Besides, this is 
good motive to scan the limit transition from smooth in- 
teraction (between BP and gas atoms) to singular "hard- 
sphere" one dividing into momentary "collisions" . 

Another our purpose is to perceive falsity of popu- 
lar treatments of the "hard-sphere BBGKY hierarchy of 
equations" (23l - [3^ trying to reduce it, in the Boltzmann- 
Grad limit, to so-called "Boltzmann hierarchy". 

We will begin in Sec. 2 by formulation of the BBGKY 
hierarchy for our particular problem and corresponding 
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(above mentioned) exact virial relations [l|, 0, [U [2^ . 
After their consideration in Sec. 3 we will see for our- 
selves that collisions constantly give birth to various sta- 
tistical correlations, and not only post-collision ones but 
also pre-collision correlations. Then, in Sec. 4, discuss 
and execute the hard-sphere limit of our BBGKY hierar- 
chy, taking in mind that it should stay compatible with 
the virial relations as they are independent on charac- 
ter of interactions (interaction potentials). The result 
is (i) modified BBGKY equations combined with (ii) fa- 
miliar boundary conditions to them at hyper-surfaces (in 
many-particle phase spaces) corresponding to collisions. 
The second ingredient allows to transform the first into 
usual "hard-sphere BBGKY hierarchy" but, at the same 
time, it forbids its reduction to the "Boltzmann hierar- 
chy" i2^|25.,ia,29.,i3i] under the Bohzmann-Grad hmit. 
Thus, the correct theory does not present such marvel- 
lous simplifications as "propagation of chaos" and closed 
equation for one-particle distribution function (DF): as 
before, one has to solve an infinite hierarchy of equations! 
This situation will be discussed in Sec. 5. 



II. BBGKY HIERARCHY, ITS CUMULANT 
REPRESENTATION, AND VIRIAL RELATIONS 

We want to consider random walk R(i) of a Brown- 
ian particle (BP) in thcrmodynamically equilibrium ideal 
gas, assuming that at initial time moment t = BP 
starts from certainly known position R(0) ~ . The 
BBGKY equations for this problem can be either de- 
rived [l^l directly from the Liouville equation, following 
Bogolyubov [H, or extracted [2l[ from general results of 
[H-i]. They reads as 



dt 



= [Hk,Fk]+n — / <i>'(R-rfe+i)Ffe+i , (1) 
ot^ Jk+i 



where A: = 0, 1, ... , is Hamiltonian of subsystem " k 
atoms -I- BP", [...,...] means the Poisson brackets, $(r) 
is potential of (short-range repulsive) interaction between 
BP and atoms, <I>'(r) = V<i>(r) , n is gas density (mean 
concentration of atoms), Fq = Fo(i, R, P ; n ) is normal- 
ized DF of the bp's coordinate R and momentum P, 
= Fk {t, R, r^'') , P^*"' in) are [k + l)-particle DFs 
for BP and k atoms [3^ , rj and pj denote coordinates 
and momenta of atoms, respectively, r'^'^-' = {ri...rfc}, 
p(fc) = {pi...p^}, and /j.... ^ J J ■■■ dvkdpk- 

Initial conditions to these equations, corresponding to 
the gas equilibrium, at temperature T , are 



Ffe|t=o ='5(R) exp(-i7fe/T) = 
5(R) Gm(P) n ■=! E{v, - R) G,„(Pj) , 



(2) 



where M and m are masses of BP and atoms, re- 
spectively, Gm{p) = (27rrr7i)"3/2gxp (-p2/2Tm) is the 
Maxwell momentum distribution of a particle with mass 
TO, and E{y) — exp [— (f>(r)/r] . Besides, existence of 



the thermodynamical limit presumes the "cluster prop- 
erty" of DFs. that is vanishing of inter-particle correla- 
tions under large spatial separation of particles, so that 
Fk — > Fk-i Gm{Ps) at jr^— R| oo , where 1 < s < A; 
and -Ffe^i docs not include and Ps . That are bound- 
ary conditions to EqlT] 

In view of these conditions, in order to visually ex- 
tract inter-particle correlations, it is convenient to make 
the linear change of variables, from the DFs Fk to new 
functions Cfc , [3J| as follow: 

Fo(t,R,P;n) = Co(t,R,P;n) , 



Fi(t,R,ri,P,pi 
= Co(t,R,P;n)/(ri~ 
+ C7i(t,R, ri,P,pi 



n) = 

-R,pi) 

n) , 



(3) 



F2(i,R,r(2),P,p(2);n) = 
= Co{t, R, P; n) /(n - R, pi) /(r2 - R, ps) + 
+ Ci(t, R, ri, P, pi; n) /(ra - R, P2) + 
+ Ci{t, R, r2, P, P2; n) /(n - R, pi) -f 

-t-C2(i,R,r(2),P,p(2);„) , 

and so on. with /(r, p) = £'(r) Gm{p) ■ 

Clearly, such defined Ck can be named cumulant func- 
tions (CF) since represent irreducible correlations be- 
tween k gas atoms and total BP's path R. In their 
terms the BBGKY hierarchy (P takes the form [ll,[l[23] 



dCk 
dt 



oP Jk+l 

Fkj Ck-i , (4) 



P d 

'mt ^ ap 



where E'{r) ~'^E{y) , pk = r*;— R, and Pfcj means 
replacement of pair of arguments Xj — {r_,-,pj} (if it is 
present) by Xk — {v}~,Pk} ■ The mentioned initial and 
boundary conditions simplify to 

Co(0,R,P; n) =J(R)Ga/(P) , 
Gfe(0,R,rW,P,pW;n) = , (5) 
Gfe(t,R,r('=),P,p('=);n) ^ at - R ^ oo 

{ \ < i < k ) . A careful enough scanning of these equa- 
tions results in observation 1221 that exact relations 



|- Gfc(t,R,rW,P,pW;n) 



(6) 



/ C,+i(t,R,r('=+i),P,p('=+i);n) 

Jk+l 



take place. This isparticular case of general "virial re- 
lations" found in [J |^, |2l| as exact properties of solu- 
tions to BBGKY equations describing molecular Brown- 
ian motion in fluids (they can be also deduced [3| from 
the generalized fluctuation-dissipation relations (ssllSGj). 

These relations demonstrate existence and significance 
of all many-particle correlations. To see straight away 
that they keep significant in the Boltzmann-Grad limit 
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(BGL), let us introduce characteristic interaction radius 
tq of the potential $(/o) , corresponding free-path length 
of BP, A — {'Kr\n)'^^ , and integrated correlations 



iyfe(t,R;A) = ^...^ Cfc(t, R, , P, p^'^) ; n) 



dV 



Thus, Wo(t,R-;A) = /Fo(t,R,P;l/7rrgA)dP is proba- 
bility distribution of the BP's path. Then Eqs|6] yield 



W^fe(i,R; A) 



d 



1 



AHaA-i 



Wo(i,R; A) 



(7) 



These exact relations hold at any values of tq and n , 
hence, under BGL (rp— )-0, n— >-oo, A= const) 
too. They show that anyway by order of magnitude 
l^fe(i,R; A) ^ Wo(t,R;A). 



III. PRE-COLLISION CORRELATIONS AND 
FAILURE OF MOLECULAR CHAOS 

Is the Boltzmann's molecular chaos ( "StoBahlansatz" ) 
compatible with virial relations ((H]), (O? May be, cor- 
relations there are purely post-collision and therefore do 
not contradict Boltzmann s ansatz proclaiming absence 
of pre-collision correlations? 

Unfortunately, this is vain hope, and the answer is no. 
To make sure of this, let us agree that post-collision { out- 
) and pre-collision (in-') states of BP and an atom satisfy 
(p-u) > or (p-u) < , respectively, - where p = r — R , 
u = V — V , with V = P /M and Vj = Pj / w being 
velocities, - and consider the first three of equations (jlj. 



9C( 



d 



dt 



9R 



dV J I 



(8) 



= — V-^^::^ + Li d — n 



dt dR 



dP 



/$'(P2)C2 -t- 

^3 



T G,n{Pl) E\pi) 



d 



dC' 



MT dV 

d 



BC 



■T(l+P21)G,„(P2)^'(P2) 



MT dV 



Co , (9) 
Ci (10) 



d 



Here, we introduced the Liouville operator 



(V 



d 

9pj 



d 



d_ 



and made use of the relative atoms' coordinates pj . 

Evidently, the last term in Eq|9] represents generation 
of pair correlation by BP-atom interaction, and the result 
is "post-collision correlation" since it passes to finishing 
"owi-state" of the particles. If we neglected three-particle 
correlation represented by CP C2 then solution to Eql9] 
would be just this post-collision correlation: Ci = Cj'"* . 
Substitution of this Ci to Eq|S] would yield a closed 
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FIG. 1. A simplest diagram transforming post-collision pair 
correlation, C""* , into pre-collision pair correlation, CJ" , 
through three-particle correlation, C2 ■ The "Brownian par- 
ticle" is represented by circles and entire arrows while gas 
atoms by dots, light lines and short arrows. 



"Boltzmann-Lorentz" equation for the BP's distribution 
Co = ^0 I thus realizing the dream of "molecular chaos" . 

However, the last term in Eq llOl for C'2 , quite simi- 
larly (even under neglect of C3 ) produces three-particle 
correlations out of the two-particle one, Cf"* . Resulting 
C2 contributes, via integral in EqlQl back to evolution of 
Ci , and produces, in particular, pre-collision pair cor- 
relations between particles going to meet one another. 
In such way, C'l acquires pre-collision component C{" , 
thus turning the dream into a rough approximation. This 
is illustrated schematically by Fig.l. 



IV. HARD-SPHERE LIMIT 

For further, first notice that EqsHllMTUlcan be written 
in equivalent form 



dCk 
dt 



V 



dCk 

dR 



d f 



+ ^ L, [Cfc + G„(p,) E{pj)Pkj Ck-i] (11) 
i=i 

Since these equations, as well as virial relations (|6]) and 
([7|) , are valid at any smooth interaction potential, we can 
extend them to the limit of "hard-sphere" interaction, for 
instance, defined by $(r) = (ro/|r|)P $0 with p — >■ 00 
and $0 — const ^ T . At that, clearly, the BP's free path 
length must tend to a constant, A ^ const . Therefore at 
any stage of this limit transition all CPs Ck are equally 
smooth functions of time (except the very beginning of 
evolution), momenta and coordinates, excluding "colli- 
sion regions" \pj\ < rQ+e , where e ^ along with char- 
acteristic thickness h of the potential wall S — ro/p ^ 
(at that, satisfying e/6 -^00). 
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In these collision regions, in opposite, Ck ( fc > ) 
become more and more sharp functions of pj , so that 
dCk/dpj ~ similarly to - 6-'^'^{pj) ■ 

Hence, the second and the fourth right-hand terms in 
Eq.(|9l), Eq. (fTO|) . etc., both are infinitely increasing (x 5^^ 
in collision regions and, in view of absence of other such 
terms, should compensate one another. Equivalently, 
the expressions in square brackets in Eqsllll should be- 
come eigenfunctions of operators "Lj corresponding to 
zero eigenvalues, that is 



L, [Cfc + G™(pj)£;(p,)PfcjCfc-i] = 



(12) 



asymptotically, inside the collision regions. 

This statement, in its turn, means that expression in 
the square brackets represents integral of motion, and its 
values at beginning and at end of the pair collision, are 
equal. Thus, visualizing only variables what change dur- 
ing collision, we can write 



Cfc(P™,p;") +G„,(p™)Pfc,Cfc_i(P™) = (13) 

= Cfc(p°"*,p;"*) +G™(p;"*)Pfc,Gfc_i(p°"*) , 

where pj = roil , with $7 being unit normal vector, and 
m -states and oMt-states correspond to J7 • Uj < and 
• Uj > , respectively, and are connected via the limit 
collision dynamics (mirror reflection): 



■r> out I out 

-t" Pj 

^ in 



(14) 



n ■ (v ™' - V ™') = - n- (vj" - V 
( 1 - r2 rj) (v/"* - V ™*) = {i-n®VL) (vj" - V ■'") 

On the other hand, integration of (IT^ over the collision 
region, - for instance, at j = k , - yields 



d_ 



<i>'{Pk)Ck 



rCk 



(15) 



((vfc - V) • 17) Ckipk - ron) dpk dn 



Substituting this equality into Eas llli for complement 
of the now forbidden regions \pj\ < vq we obtain 



dCk 
dt 



-V-^-J2ur^ -nTCk+, , (16) 



dR 



dpj 



where Uj = Vj — V . At boundaries of the forbidden 
regions, i.e. at \pj\ — ro , these equations should be sup- 
plemented with boundary conditions p^ - ([T^ . Combin- 
ing Eqs[T6]with ((T5|) and ([13]) we come to 



ttA 



dCk 
dt 



dp 



OR 



dfl ((V - v) • SI) X (17) 



(V-v)-0> 



X [G„(p*)Cfe(P*)- G„(p)Gfe(P) + 
- Ck+iip = ron,P*,p*)- Ck+iip = roQ , P, p) ] 



with p = pk+i, p = Pfc+i . V = Vk+i = Pfe+i/m, 
and momenta P* , p* being related to P , p in the 
same way as pre-coUision momenta P , p ™ in ([H]) 
are related to post-collision ones, P°"* , p°"* , and the 
integration involves pre-collision states only. 

What is for the operators L° , in respect to CFs they 
are defined as L° = — (vj — V) • d/dpj at \pj\ > ro 
and by boundary conditions (|13p at \pj\ — tq . Thus, 
importantly, are not mere translation operators: at 
collision surfaces \pj\ = tq they act as creation oper- 
ators, creating (k + l)-order correlations from fc-order 
ones. Therefore factually Ck remain coupled with both 
Gfe+i and Ck-i , as in basic EqsU) In other words, 
due to conditions p^ . hierarchy ([T7| keeps character- 
istic three-diagonal structure of equations for CFs ! 

Return from CFs back to DFs, according to the CFs 
definition ([3]), transforms these equations to 



dFk 

at 



dFi 



dFk 



H-r^-nVFk,, (18) 



and conditions (IT^ to 



Ffc(P"\pj") = Ffc(P°"*,p;"*) at |p,|=ro (19) 

along with Their substitution into the "collision 

integral' T Fk+i gives, similarly to (IT71) . 



ttA 



dFk 
dt 



dp 



dR 



dVL ((V X 



(20) 



(V-v)-Ji>0 



X [Fk+i{p = rnn,l?\p*)-Fk+i{p = ron,V,p)] 

Here again = — (vj —V) • d/dpj at \pj\ > tq but at 
\pj\ — tq action of the operator L° onto DFs is defined 
by boundary conditions ([1 



THE BOLTZMANN-GRAD LIMIT AND 
"MATHEMATICAL NON-PHYSICS" 



Equations (|16p as combined with boundary conditions 
P^ - fTI)) at collision surfaces \pj\ — tq (and conditions 
([5]) at infinity) or, equivalently, equations (1181) combined 
with (fT9|) and (fT4|) represent direct analogue of so-called 
"hard-sphere BBGKY hierarchy" (see e.g. Eq.2.2 from 
[H] or Eq.4.1 from '^). 

Importantly, the term "hard-sphere BBGKY hierar- 
chy" is adequate on the understanding only that when 
the boundary conditions are substituted into Fk+i inside 
the "collision integral" in equation for Fk then simul- 
taneously and necessarily they are satisfied by Fk+i 
in the next equation for Fk+i itself (i.e. included into 
definition of the operators L° as above) . Otherwise one 
makes some "mathematical non-physics" , since resulting 
equations will be non-derivable from Liouville equation 
for an (infinitely) many hard sphere system ! 
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In opposite, fulfilment of the mentioned requirement 
guarantees observance of the virial relations (see Ap- 
pendix) and consequently belonging of resulting equa- 
tions to the class of BBGKY hierarchies (since virial re- 
lations do follow already from most general properties of 
many-particle dynamics and Liouville equations [H-Q)- 

In view of these facts, it seems impossible to accept 
the old idea that under the Boltzmann-Grad limit (BGL) 
the "hard-sphere BBGKY hierarchy" is equivalent to so- 
called "hard-sphere Boltzmann hierarchy" (see e.g. [l^, 
HI) m, mi). For our system it looks as 

dt dR ^ ' dpj 



pre-collision correlations do arise, at \ij -pj < and again 
at |uj X pj|/|uj| < rg , since statistical picture acquires 
more and more time-reversal symmetry. 

At the same time, at |uj x pj\/\uj \ > rg there are al- 
most no reasons for correlations. Hence, in directions 
perpendicular to Uj all CFs become more and more 
sharp functions of pj . At that, the integral value of 
correlation per one atom goes to zero oc TrrgA = 1/n. 
But this has no significance since correlations concentrate 
just where they are most effective. This is confirmed by 
the virial relations ([T]). The latter show also that really 
significant correlational characteristics are n^Ck which 
turn under BGL into singular generalized functions. 



9R 



(V-v)-0> 



dn ((V- v) • n) X (21) 



VI. CONCLUSION 

We have seen that the Boltzmann-Grad limit does not 
lead to any essential simplifications (and in this sense it 
does not exist). In particular, it does not lead to the 
mythologic Botzmannian kinetics with Boltzmann equa- 
tion, Boltzmann-Lorentz equation, etc.), except non- 
interesting spatially homogeneous case [39[ . Physical and 
statistical reasons of this pleasant fact already were over 
and over again explained in [^-[1] and earlier |T2l - [l7| and 
later P, i, i, [H [II EH . The fact is pleasant because 
it shows once again that real many-particle dynamical 
chaos does not degenerates into miserable stochastics. 



Integration of (fc+ll-th of EqslTBlover Pfc+i and Pk+i 
at |pa;+i| > '"0 yields for Ck = /^.^i Ck+i , in view of (O, 
equations 



PCfc+i - nPCfc+i 



x[i^,+i(p = 0,P*,p*)-Ffe+i(p = 0,P,p)] , 

where there are no forbidden regions, that is pj can take 
arbitrary values from the whole R'^ , and the boundary 
conditions (IT5|) at |pj| — J> rg ^> are thrown, that is op- 
erators are replaced by trivial translation operators, 

-(v, -v). a/ap,.^ 

This means that in corresponding equations for CFs 
any Ck is now connected to Ck+i only, and not to 
Ck-i ■ Thus, these equations form two-diagonal hier- 
archy qualitatively different from the original one ! 

As the result of such frivolity, these equations allow 
for factored solution with "propagation of chaos" , when 

Fk = Fq 11^ = 1 G'„i(pj ) , Cfc = at fc > . and Fq — Co Appendix: Virial relations for hard-sphere system 

undergoes the "Boltzmann-Lorentz equation" . 

Evidently, this "Boltzmann hierarchy" contraries to 
the above emphasized requirement: it uses the boundary 
conditions to write "collision integrals" but neglects the 
same conditions in higher equations what determine the 
integrands (as if the latter took particles from a "parallel 
world" ) ! As the consequence, the Boltzmann hierarchy 
does not satisfy the virial relations. Its wrong is clear 
already from observation that it results when one first 
damages Eqs HO] by replacing with —{\-j—'V)-d/dpj 
and only after that goes to ro — 37]. Therefore, it is 
not surprising that the famous Lanford's attempt |24l to 
prove the mentioned idea was in fact unsuccessful 8, 38]. 

Possible formal cause of this nonsuccess is very simple: 
the smaller is rg the less smooth functions of pj are all 
the CFs since conditions ensure continuity and 

smoothness of the probability measures Fk along phase 
trajectories only but not in perpendicular dimensions. 

Indeed, let Uj ■ pj > , that is BP and an atom scat- 
ter one from another. If at that rg \pj\ < A and 
additionally \uj x pj\/\uj\ < tq (where x denotes vec- 
tor product), i.e. vectors and pj are nearly paral- 
lel, then the particles form an out-state which (almost 
certainly) arose from their recent collision (especially if 
\pj\ ^ A). According to (IT^ . post-coUision correlations 
of any order between such particles inevitably take place, 
and by order of magnitude all they are equal to Cg = Fg . 
With time, as was explained in Sec. 3, also quite similar 



dCk 
dt 



dCk 



i=i 



dCk 
dpj 



At the same time, differentiation of Eas ll6l in respect to 
the density produces for Ck — dCk/dn equations 



dCk 
dt 



dCk 
' dK 



k 

E 

i=i 



dCk 
dpj 



PCfc+i — nVCk+i 



which are identical to the previous ones. Since, according 
to (UJ, initial conditions for Ck and Ck also are identical 
(all are zeros), we come to equalities Ck — Ck ^ that is 
to the virial relations dU, dCk/dn = Ck+i ■ 

Notice that this consideration is valid regardless of con- 
crete form of boundary conditions at \pj\ = rg . There- 
fore, virial relations hold under replacement of the hard- 
sphere "mirror reflection" conditions ([H]) by some other 
rule. Correspondingly, equations ()16p and ()18|) can be ex- 
tended to collision operators P with other "scattering 
laws" than (IE 
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